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It is pointed out that electron thermal fluctuations can couple with the ion 
acoustic mode in an inhomogeneous plasma to generate a low frequency ion time 
scale electromagnetic wave. This electromagnetic wave can become unstable if 
the temperature and density gradients are parallel to each other which can be 
the case in laser-plasmas similar to stellar cores. The comparisons of the present 
theoretical model with the previous investigations are also presented. The final 
result is applied to a classical laser induced plasma for illustration. 

I. Introduction 

A low frequency electromagnetic wave should be a normal mode of unmagne- 
tized inhomogeneous plasmas as an intrinsic source of magnetic fields. The ion 
acoustic wave (lAW) is a very important mode because it exits in both unmag- 
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netized and magnetized plasmas. It is a low frequency wave and light electrons 
are generally believed to act as inertia-less in its electrostatic field. 
In fact the electron inertia can play a vital role in producing a low frequency 
electromagnetic wave. More than a decade ago [1] , it was proposed that an elec- 
tromagnetic wave having frequency near lAW is a normal mode of unmagnetized 
plasmas which can be responsible for magnetic field generation if it becomes un- 
stable. In the derivation of its linear dispersion relation the electron inertia was 
not ignored while the displacement current was neglected. It was shown that 
the compressibility and vorticity can couple due to density inhomogeneity and 
hence a low frequency electromagnetic wave can be produced. Basically elec- 
trostatic lAW and magnetostatic mode [2] cooperate with each other to develop 
such a wave. The electron temperature perturbation was not taken into account 
and the steady state was assumed to be maintained by external mechanisms. 
The theory was applied to explain the magnetic field generation in laser plas- 
mas. But the longitudinal and transverse characters of electric field decouple if 
the quasi-neutrality is used. Then both the high frequency and low frequency 
electromagnetic instabilities were also investigated in unmagnetized plasmas [3] . 
The interest in the investigation of low frequency magnetic fluctuations in un- 
magnetized plasmas was initiated after the first experimental observation of 
magnetic field generated in a laser induced plasma [4]. Then more experiments 
were performed [5, 6] on these lines. Several theoretical models were presented 
to explain the magnetic field generation in laser produced plasmas [7-13]. 
The magnetic electron drift vortex (MEDV) mode was proposed as a pure trans- 
verse linear mode which can exist because of the electron temperature fluctu- 
ations in immagnetized inhomogeneous electron plasmas [13]. This mode can 
become unstable [14] if the equilibrium electron temperature gradient is parallel 
to the density gradient and it is maintained by external effects. 
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Most of the mechanisms proposed to explain magnetic fluctuations in initially 
unmagnetized plasmas are based on electron magnetohydrodynamics (EMHD) 
which has been discussed in detail in Refs. [15, 16]. 

Recently [17], EMHD has been used to study long-lived and slowly propagating 
nonlinear whistler structures (NLWS), or whistler spheromaks (WSPS) observed 
in laboratory experiments [18, 19]. In AIEDV mode the divergence of electric 
field is assumed to be zero (V.Ei) = while the divergence of electron velocity 
is non-zero (V.Vgi ^ 0). Furthermore the electrons are inertia- less and ions 
are treated to be stationary. The frequency w of the MEDV mode is assumed 
to lie in between the ion plasma frequency and electron plasma frequency i.e. 

<< a; << where = , for j=e,i and c is speed of light, while 

k is the wave vector. These restrictions and assumptions are indeed very strict. 

It is important to find out some electromagnetic mode taking into account the 
ion dynamics so that the strict restrictions on the frequency and wavelength of 
the perturbation are relaxed. Then the only required condition on frequency 
should be w << ujpg.ck. Furthermore, the electron temperature gradient is 
necessary for hydrostatic steady state VpeO = 0. In this case the temperature 
gradient becomes anti-parallel to density. But laser and astrophysical plasmas 
are open systems and many external mechanisms can maintain a study state 
with parallel density and temperature gradients. For example, in stellar cores, 
both the density and temperature increase towards centre of the star due to fu- 
sion and star is held intact because of gravity. Therefore both the parallel and 
anti-parallel gradients can be discussed. The electron thermal fluctuations can 
produce electromagnetic wave as was proposed many decades ago [13]. Several 
interesting works have appeared in literature on these lines (see for example 
review articles [9, 15, 16]). 
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It is pointed out here that it is not necessary to assume a pure transverse 
perturbation in an electron plasma. Rather the perturbation can be partially 
transverse and partially longitudinal and hence the electron density perturba- 
tion may not be neglected. The frequency of such a wave in electron plasma 
turns out to be near VteKn where v^, = (^y^^)^ if' the electron thermal speed 
and K„ = the inverse of density gradient scale length L,,, = In 

the local approximation we need to have k„ << k. Furthermore the condi- 
tion << ulgfc^ can be satisfied if ^ << X%gky. Therefore, generally we 
may have vteKn < ojpi and hence one cannot neglect ion dynamics. Moreover, 

k (where Cs = ' is ion sound speed) is also possible. 
Therefore, it is expected that the electron thermal fluctuations can couple with 
lAW to produce stable and unstable low frequency electromagnetic waves in 
unmagnetized plasmas. We are interested mainly in finding out a quadratic 
dispersion relation for a low frequency electromagnetic wave to see a simplest 
physical picture. For a comparison and general interest, the previous theroreti- 
cal models are also briefly discussed. 

In the next section, the pure electron plasma is investigated including the de- 
scription of MEDV mode. In section III, the coupling of lAW with magnetostatic 
mode due to inhomogeneity is briefly revised. In section IV a new electromag- 
netic mode which arises due to the coupling of lAW with thermal fluctuations 
is presented.In section V the result is applied to a case of classical plasma for 
illustration. Finally the physical assumptions and mathematical limits of the 
results are discussed. 

II. Electron Plasma 
Let us consider the electron plasma in the background of stationary ions. 
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The set of equations in the linear hmit can be written as, 



memdtVei = -enoEi - Vpei (1) 



V X Bi = — Ji (2) 

c 



Ji = -enoVel (3) 



V X El = --^tBi (4) 

c 



Since 



pi = noTei 
therefore energy equation becomes, 



3 3 

-nodtTei + -no(Vel.V)TeO = -poV.Vel (5) 



First we describe the pure transverse MEDV mode. 
Curl of (1) gives. 



at(Vxvei) = — atBi + ^— (Vno)xVTei (6) 

nieC meUo 



Equations (2) and (3) yield, 



and hence 



Vei = -T-^V X Bi (7) 
47reno 



V X Vel = V'Bi (8) 



where Vno x (V x Bi) = due to the assumption k_LVno-LBi. Equation (8) 

dnp 
dx 



predicts Ei = Eix while Vno = V = {0,iky,0) and Bi = Biz have been 
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chosen. 

Equation (6) and (8) yield, 



(1 + A^fc2)atBi = (Vno X VTei) (9) 

eno 



where Ag = Equation (7) gives, 



V.Vei = -^^.(VxBi) (10) 

AimQC no 



and therefore one obtains. 



2 c 

Tel = ^-j kyKnBi (11) 

3 i-KTlQC 



Then (9) and (11) give the Unear dispersion relation of MEDV mode as. 



= \Co{^fvhkl (12) 



where Cq = i+'xik'^ ^'^^ = {Tf./me)2. If VTeo ^ is assumed, then (12) 
becomes, 



2 r ^-n 



ky 



where kt = 1 7^^^ I and VTeo = has been used. If (|k„ - kt) < 0, 

then the mode becomes unstable [14]. 

It can be noticed that we have V.Ei = and V.Vei 7^ in the above treatment 
and it does not seem to be very convincing. 

If we write (1) in x and y components, we can observe that the y-component of 
El should not be considered as zero due to dyPei ^ and hence the wave may 
be partially transverse and partially longitudinal. 
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The curl of (2) yields a relation between Ei^; and Eij^, 



(14) 



where a = (1 + A^fc^). 
Equation (5) becomes, 



W, 



2±el 

To 



2kl 



nei 
no 



e 

Tor 



ikyEiy + -krEix 



(15) 



where Wq 



3/ j2 

5'^ 



vlkl ( 1 



(l-i-^) 



The continuity equation yields, 



1 + I I i kyEly 



me 



-— 1 + ^ 



TOe \ 2 Kn 



r KnEu 



where = jw^ - v^^k^ - (l - f^) }• The Poisson equation 



(16) 



V.Ei = -47re 



"-el 

no 



(17) 



can be written as, 



vlkl(^ + Gn)u^'-^vikt{g'T-9lT) 



pe 



3 2 I 2 7 2 2 
0^ +7^'"teky9nT 



iE- 



■"2 2 I 2 



-,9nT - 



3 kt 

2 Kyr 



(18) 



where G„ = (l - §,9^^), = and gf, = Kf^/fc^. 

Equations (14) and (18) yield a linear dispersion relation in the limit << 
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as, 

(19) 



2 2 

OJ — 



where i^o = {(1 + A|,,fc2)a - f| + a\\^klG^^ and A|,, = (4;^) • The last 
term on right hand side of (19) is smaller because Afj^fc^ < 1 and due to local 
approximation we have k„, kt « ky. Therefore (19) can be simplified as, 



which predicts a purely growing electromagnetic mode. 

The dispersion relation (13) can not be retrieved from (20) because we have 
evaluation V.vie from equation of motion and in the derivation of (13), the 
velocity Vei is ovahiatcd using Ampere's law for pure transverse wave. The 
frequency uj in both the dispersion relations is of the order of vf^K^ or smaller. 
The important point to note is that 

and Xjjgky « 1 has been assumed to ignore electron plasma wave. Therefore, 
even for ^ < ^f^^ , the frequency in (21) can satisfy the relation to < Wpi or 
even w can be of the order of lAW frequency Cgky. 

Therefore, the partially transverse and partially longitudinal electron wave of 
equation (20) may couple with ion acoustic mode. 



III. lAW and Magnetostatic Mode 

Here we shall show that transverse magnetostatic mode [2] which is obtained in 
the limit uj^ << lo^^ can couple with lAW in a nonuniform plasma [1]. Therefore, 
both ions and electrons are considered to be dynamic. The ions are assumed to 
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be cold for simplicity and therefore equation of motion becomes, 



9tVii = —El (22) 



The continuity equation yield, 



Tl"\ €■ 

— = 2 ^i^nExa^ + i^Eiy) (23) 

no rriiU)'^ 



If quasi-neutrality is used then transverse component Ei^ and longitudinal com- 
ponent E\y become uncoupled. Therefore, we assume Af^gA;^ ^ and use Pois- 
son equation which in the limit << becomes, 

(24) 

Note that the term vj^ky is not ignored compared to uj^^ to couple Ei^ and Eiy. 
Equations (14) and (24) give a linear dispersion relation as. 



.2 



ui 



{ah - ^) 



(25) 



where = and 6 = (l + Alj^fc^). We have to use Poisson equation to obtain 
a quadratic equation in uj while (2) implies nei ~ n^i. The equation (25) is the 
same as equation (21) of Ref. [1] where two small terms in the denominator are 
missing. Note that A|,g < Ag and if quasi neutrality is used due to Ampere's 
law, then (25) yields the basic electrostatic lAW dispersion relation = (?gky. 
If displacement current is retained and Poisson equation is used without using 
« Wpg, then one obtains a dispersion relation of coupled three waves; 
ion acoustic wave, electron plasma wave and high frequency transverse wave [3] . 
Actually the contribution of displacement current in the curl of Maxwell's equa- 
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tion has been neglected for << uj^^^c^k'^. This should not mean that the 
electrostatic part of current is also divergence free, in our opinion. In the diver- 
gence part of Maxwell's equation << is used but v^^ky term is assumed 
to be important. It may be mentioned that in this treatment, Ampere's law 
does not imply quasi-neutrality necessarily. We need a coupling of divergence 
part and cure part of the current and for this we need to assume 1 < X^^ky in 
the limit 1 < << Wpg,c^P. 

IV. Unstable Electromagnetic Wave 

Now we present a simple but interesting picture of ideal plasma assuming ions 
to be cold. The divergence and curl of (1) give, respectively, 

atV.(noVei) = — ^noV.Ei - — Vno.Ei - — (V.Vpei) (26) 

nie TOe me 

and 

9t(V X Vel) + (Kn X dtVel) = —Kn X Ei - —V X Ei (27) 

rUe me 

where k„ = \^^-\ and Vno = has been assumed. If initially electric 

field was purely electrostatic i.e. Ei = —Vipi, then it will develop a rotating 
part as well if Vno x Ei ^ 0, as is indicated by the right hand side (RHS) of 
(27). 

Therefore, in an inhomogeneous plasma electron thermal fluctuations can take 

place due to low frequency perturbations of electric field which has both trans- 
verse and longitudinal components. Then these thermal fluctuations can couple 
with I AW even in the quasi-neutrality limit. 
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However, we use the Poisson equation 



V.Ei = Annoe 

V no no J 



Using (16) and (23), the above equation can be written as. 



In the hmit << tOpg, (29) reduces to. 



rSo 4 \ 1 1 1 1 I ^ '^2 272 2 2 7 2.^ i 



Equations (14) and (30) yield a linear dispersion relation as follows. 



+ (^) ' [I + G„) + vlklGr. - Ivlk^!^]].^ 
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+laclklvlkl{9'^ - 9l^) - [^fj clklvlklial - gl^) = (31) 

This equation contains a coupling of both the electromagnetic waves described 
by the dispersion relations (20) and (25). We are interested in finding out a 
simple dispersion relation in quadratic form to look at the physical phenomena 
more clearly. 

For this we assume g„T = gr i-e. kt = and in this case (31) becomes, 



3Ho 



[cX + G.) (a - Kl/kl) - vUl (l + Ixlkij I (32) 



where Hq = 



In the quasi- neutrality limit, we have 



b=l in the above electromagnetic dispersion relation and it reduces to 

where iJi = (o — K^/ky) and for «;„ = 0, we obtain = ^c^ky. 
It is important to mention that the transverse and longitudinal components of 
electric field do not decouple in the presence of electron thermal fluctuations 
even if quasi-neutrality is used. If thermal fluctuations are ignored (32) reduces 
to (25) and if ion dynamics is neglected we recover the result of equation (20). 
In principle both the waves described by (20) and (25) have frequencies very 
near to each other under local approximation in plasmas of low-Z materials. 
Therefore the coupled mode (32) is a fundainciital unstable electromagnetic 
wave of unmagnctizcd inhomogeneous plasmas. If X^ky < rrie/mi, then a ~ 1 
and b=l. In this case Eq. (33) predicts purely growing electrostatic fluctuations 
for ^ < (f-) because the second term on RHS of (33) can become larger 
than the first term. 
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V. Application 

As an illustration we consider the parameters of classical laser-induced plasma 
as no ~ 10^°cm~^ and Tg = lOOev. Then we obtain Cg ~ 10^ cm/ Sec, Vte = 
4.18 X lO^cm/Sec, Ae ~ 0.53 x IQ-'^cm, Xoe ~ 0.74 x 10-^cm, and u)pi ~ 
1.3 X lO^'^rad/S'ec. If k„ ^ lO^cm^^ and wc choose ky ^ lO'^cm^^, in Hy- 
drogen plasma wc find ^ < (^^Y = IQ-^ . Wc obtain a=29, b=1.0055, 
Vtei^n ~ 4.18 X lO^'^rad/Sec and Csky ^ lO^^rad/Sec. Since ^ < Xj^^ky, we 
use the result of Poisson equation (32) and obtain 



Note that it gives a purely growing electromagnetic instability. Furthermore 
Cgky < \u}\ < u)pi and hence one should not ignore ion dynamics. 
The instability appears due to the second term on RHS which enters because 
of pure electron dynamics. But fluctuations are at ion time scales. Therefore, 
ions must be treated to be mobile. The instability can also appear in the quasi- 
neutrality limit. 

If Kn ^ 10"^ cm^^ and we choose ky ^ lO^cm"^, then X^^ky < nie/rrii and hence 
we use quasi-neutrality. Then a=1.28 and (33) yields, 



Again we see a purely growing electromagnetic mode with < u>pi. 
VI. Discussion 

A low frequency unstable electromagnetic wave has been studied in unmagne- 
tized inhomogeneous plasmas. This wave can be an intrinsic source of magnetic 



oj ~ ±i(5.8 X 10^2) 



(34) 



Lu ~ ±i{3.6 X 10^2) 



(35) 
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field fluctuations in initially unmagnetized plasmas. Dynamics of both electrons 
and ions are important in this instability. 

These fluctuations can become unstable when the density and temperature gra- 
dients are parallel to each other as is the case, in the stellar cores. Similarly 
in laser-induced plasmas these gradients can be in the same direction due to 

external conditions. 

In a statistical ensemble of charged particles, the electrostatic fluctuations arc 
quite natural. Any initial electrostatic field perturbation can produce its trans- 
verse component as well in the presence of density gradient. This phenomenon 
can cause a coupling of ion acoustic wave (lAW) with the low frequency trans- 
verse magnetostatic mode. This coupled mode has already been investigated 
more than a decade ago [1]. But it can exist in a relatively shorter wavelength 
range for ^ < Xjj^k'^. In the quasi- neutrality limit, the I AW and magneto- 
static modes decouple. 

Another interesting mechanism has been proposed for the generation of a trans- 
verse wave in as electron plasma. In this case the transverse wave (MEDV mode) 
is produced by the electron temperature perturbation [13, 14]. This mode does 
not have density fluctuations. In the present investigations it has been shown 
that a low frequency electromagnetic wave seems to appear in a non-uniform 
unmagnetized electron plasma having both the longitudinal and transverse elec- 
tric field components. But the frequencies of both the above mentioned modes 
are of the order of Vtei^n or even smaller. 

Since frequencies VteKn and Cgky can be very near to each other, therefore ions 
should not be assumed as stationary. We need to consider ion dynamics as 
well. Then we obtain a coupled linear dispersion relation which is fourth order 
polynomial in lu. If = k„ is assumed, then it reduces to a quadratic equation 
(32). It indicates the presence of a partially transverse and partially longitudi- 
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nal unstable wave in nonuniform unmagnetized electron-ion plasmas. Dynamics 
of both electrons and ions take part in the generation of this instability. The 
electron temperature fluctuations are necessary to produce this wave. It is also 
important to mention that this wave can exist in the quasi-neutrality limit as 
well. 

In our opinion, this electromagnetic instability plays an important role in the 
generation of fluctuating magnetic fields in initially unmagnetized plasmas. The 
nonlinear saturation mechanisms need to be investigated which can contribute 
to the magnetization of the plasmas. Therefore the electromagnetic waves and 
instability presented here can be important in the study of magnetic field gen- 
eration in laser-induced and astrophysical plasmas. 
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